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Abstract 



We illustrate some of the techniques to identify chaos signatures at the 
quantum level using as a guiding examples some systems where a particle 
is constrained to move on a radial symmetric, but non planar, surface. In 
particular, two systems are studied: the case of a cone with an arbitrary 
contour or dunce hat billiard and the rectangular billiard with an inner 
Gaussian surface. 
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1 Introduction 

Hard wall billiards are among the simplest and most studied systems in the field 
of classical and quantum chaos [I] -[5]. A two dimensional quantum billiard is 
analogous to a vibrating membrane because of the mathematical equivalence 
of the stationary Schrodinger and Helmholtz equation with the same Dirichlet 
boundary conditions. This similarity between billiards and classical waves has 
been used in experiments e.g. the quantum microwave cavities [6]. When the 
classical counterpart of these systems are chaotic, some signatures of this chaos 
appear at the quantum level. One of them is the statistical properties of the 
spectrum, which follow the Bohigas, Giannoni and Schmit conjecture [7]. This 
conjecture states that the nearest neighbour energy level statistics of a time 
reversal invariant chaotic system follows the same statistical as the ensemble 
of random orthogonal matrices with Gaussian distributed elements, i.e. the 
Gaussian Orthogonal Ensemble (GOE) [8 . For 2x2 matrices this distribution, 
also known as the Wigner distribution, is given by 



To be more precise, the energy levels should be first classified according to 
the partial symmetries that the system has, then the statistical analysis should 
performed on each symmetry class of states, and for each of those the nearest 




(1) 



neighbour spacing distribution of the energy level, properly normalized, follows 
the GOE distribution ([T]). On the other hand, for systems with enough sym- 
metries to be integrable, the nearest neighbour energy level statistics follows a 
Poisson distribution. 

Additionally, quantum chaotic billiards exhibit the phenomenon of scarring 
of their wavefunction specially when the corresponding state is deep in the semi- 
classical limit 0. Although, the scars in the wavefunction were first observed in 
the seventies [10] they have become a phenomenon of great interest because they 
represent a connection between the quantum chaotic billiard with its periodic 
orbits. 

The chaotic features of two dimensional billiards are often studied confining 
them on a plane with a contour dD where chaos emerges when an "irregular" 
contour is used. However, the system present interesting features if we set a 
non-planar geometry in the billiard internal region ^, even when the contour is 
regular as a circle, a rectangle or an ellipse. 

The aim of this work is to illustrate some of the characteristic signatures of 
chaos at the quantum level in non-planar billiards. As such, this paper is an 
introduction for graduate and advanced undergraduate students to the topic of 
chaos in quantum systems, illustrated in systems different than the usual planar 
billiards. It can serve as a guiding example to the techniques used to analyze 
and identify chaos signatures in quantum systems. Two systems are studied: 
the dunce hat billiard (or conical billiard) and the rectangular billiard with a 
Gaussian surface. These billiards are shown in figure [l] 




Figure 1: The systems: (left) the dunce hat billiard with an arbitrary contour, 
(right) the rectangular billiard with a Gaussian surface immersed in a uniform 
electric field. 

We start in section |2] with a geometrical background needed to study non- 
planar systems. Although most of this work is about quantum systems, in 
section |3) we study the classical dynamics on these billiards and present some 
of the chaos characteristic signatures at the classical level. Then we proceed. 
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in section |4j to present how the quantum Hamiltonian is constructed on non- 
planar surfaces, and in sections |5] and [6] both types of billiards are studied at 
the quantum level. 



2 Geometrical background 

Consider a surface S embedded in the three-dimensional Euclidean space Es 
defined by 

^:y' =y\u\u^), (z = l,2,3) (2) 

where {y^ ,y^} are the Cartesian coordinates in and {u^^u^} are the 
Gaussian curvilinear coordinates defined on the surface. A particle restricted 
to move on a two-dimensional surface has only two degrees of freedom, and its 
position r G may be expressed as follows 



(3) 



The motion of a particle will depend of the surface geometry which is specified 
by the first and second fundamental forms. Several of the intrinsic properties of 
H are determined by its metric aaf3 included in the first fundamental form ^/ 
given by 



:= df'df= aasdu^du^ with aas = 



Qyl Qyl 



(i = l,2,3)(a,/3 = l,2) . 

(4) 



We are interested in billiards with a radial symmetric surface defined by 

y'^ = cos {u^) y'^ = sin [u'^) y^ = f{u^) . (5) 



For this particular case the metric is a diagonal tensor 



(6) 



Therefore, the system is orthogonal and the Gaussian curvature K may be 
computed with 



K 



1 



d f 1 da 



22 



d f 1 dan 



2y/a \_du^ J du^ Vv^ 

where a := det^aap)- Replacing the metric given by the equation (|6| we find 



(7) 



K ■■ 



(8) 
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Classically, the motion of a particle which lives on S is affected by the surface 
curvature because the surface metric is included explicitly in the kinetic energy, 
which is proportional to aa^ii^u^- At the quantum level, this kinetic energy is 
associated with the Laplace-Beltrami operator on the surface. Other extrinsic 
properties of S are specified in its second fundamental form ^jj given by 

1 / dr dr \ d^r 
^n:=b^,du-du^ with b^, = — ^—^—y^-^, (9) 

and the matrix (ba^) for a radial symmetric surface is 



This enables to compute the principal curvatures ki and k2. If the principal 
directions coincide with coordinate curves, then 

ki = and k2 = . (11) 

an a22 



Finally, if the surface is radially symmetric, then 

ki = ^ ^ and k2 = jM^=. (12) 



3 Classical motion 

The Lagrangian in cylindrical coordinates of a particle of mass ji with a time 
independent potential V is 

jSf = |(r^ + ^2^2 ^ .2^ _ y^^^ (^3) 

The position of particle constrained to move on a radial symmetric surface 
defined by z := /(r) is r = r (1/1,1^2)- Then, the Lagrangian takes the form 

^=^a^pu''uP -V{r) (a,^ = l,2) (14) 



with 

= {hpf with (/ll,/l2) = ^^1 + {drff^V^ , (15) 

so the canonical momentum is 

Pa = duc^ = fiaa^u^ . (16) 
Therefore, the Hamiltonian is 

= pau^{u,p)-^ = ^a'^^PaP^ + V (f) with a^^ = {h^)''^ . (17) 
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If the particle has a charge ^ < 0, and the system is placed in a uniform electric 
field E = Eok (with Eq> then V(r) = -qEof{r). Finally, the Hamiltonian 
is 

and the Hamilton equations of motion for the r-coordinate are 

Pr 



Pr = —drJ^ = 



1 + {drfY 
I Pr 



qEo 



drf- 



3. 



(19) 



(20) 



The z-component of the angular momentum remains constant in between two 
successive collisions of the particle with the billiard walls. As a result, for the 
(7!)-coordinate 



9p^^ = — 2 



with Lz = if r{t) < rdcj)) , 



(21) 



where the billiard contour is 9!D = {(rc ((/>), (j)) - < (j) < 27r} and L^^ is a con- 
stant. The value of L^^ depends of the billiard contour. For instance, L^^ has 
the same value for all the collisions in a circular contour but in general it must 
be updated after each collision when the contour is not a circular one. 



3.1 The dunce hat biUiard map 

After the plane surface, one of the simplest case of study with zero Gaussian 
curvature is the cone. Let us define the dunce hat billiard as a cone /(r) := 
— -^) with a circular contour = {(i?, (/>) : < < 27r}. The Hamiltonian 
of this billiard in absence of the electric field according the equation (18) is 



pI 



2/i 



2/ur2 



(22) 



This Hamiltonian differs from the free particle on the plane surface Hamiltonian 
for the constant term , and avoiding the circular contour the free particle 

solution on the cone is found by direct integration of the equations of motion. 
The position is given by 



r(t) 



i 



1 



(ft)- 



(23) 
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and 



— arctan 




where (ro,^o) is the initial position of the particle in polar coordinates. The 
radial and angular velocity at t = are and uJq respectively. Although, the 
structure of r{t) and (j){t) on the cone is similar to the one obtained for the free 
particle on the plane surface, the cone may deflect the particle trajectory as it 
is shown in Figure |2] 
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Figure 2: Trajectory deflection of the free particle living on the cone. 

This plot is a top view of the trajectories followed by a particle when the cone 
height varies from = to = R. The typical rectilinear trajectory of the 
free particle on the plane /o = is deflected in a counter-clockwise sense as 
the cone emerges by increasing fo specially near the origin where the cone is 
located. 



arctan 




2 ToVr^ + 



The solution of the equations of motion on the cone, including the circular 
contour, may be found defining the following variable rescaling 



T : r{t) 



and 4){t) = \f^o(t)(t) with := 



1 



(25) 



The Lagrangian with the new variables takes the form 



{4y 



(26) 
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with X = fcos(^ and y = rsin^. The new canonical momenta are: 

Pr V^Pr and p0 = ^|=. (27) 



By computing the Poisson brackets of the new coordinates it can be checked 
that the transformation is canonical. 



Q 11 







Figure 3: Classical trajectories of the circular dunce hat billiard, (left) 
Periodic orbit on the XY frame, (center) Mapping from the XY to the XY 
frame, (right) Trajectory of the particle on the surface. We have set fo = ^ 
and fo = ^ for the upper and lower panel respectively. The trajectories in the 
circular billiard are periodic orbits. Note that cone may generate orbits which 
are not periodic. 



The region D remains as a disk: x'^ 



contour is 



= x^ ^ y'^ < . However, the 
< < ; ^ > . The upper limit of 6 may 

be written as c;^ < V^^tt = 27r — a where a := (l — V^) ^tt. If we take a 
circle of radius R' := \/ -\- R^^ then we may build a cone of radius R and 
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height fo by removing a circular section with angle 27r 



As a 



result, a is just the angle related with the circular section removed from the 
circle in order to construct the cone. In the XF-plane we have a circular planar 
billiard. Therefore, if (j)n locates the n-collision and /3n is the angle of the incident 
velocity of the n-collision, then the next collision is connected with the billiard 

map Mn : ^^n^/^n^ (^^n+i7 /^n+i^ This map for a circular billiard is given 

by in 



9n+l 



1 2 
1 



eir 




(28) 



with (e = +1 if the motion is clockwise, -1 otherwise). Using the equation (28) 
the particle position is 



x{t) = R 



(^COS^^(^) + l -COS^^(t)) (J^ -n{t)^ +COS^^(^) 



and 



y{t) = R 



(sin ^n(t)+i - sin ^^(^)^ (- - n(t) j + sin ^ri{t) 



(29) 



(30) 



Where r is the time between successive collisions, and n := Int(t/r) is the 
integer part oft/r. According to the circular billiard map = ^o+^(2/^o + e7r). 
Therefore, at the n-th collision the particle is located at 



n 




with Pri = Po ' 



(31) 



The last expression gives us the position of each collision. However, this map 
does not specify how successive collisions are connected. As we show in Figure 
|2]the cone deflects the particle. Therefore, two successive collisions predicted 
by ( [31] ) must be connected by a curve (see Figure [3|. This is an expected 
difference with the plane billiards where collisions are connected by rectilinear 
trajectories. 



3.2 Classical rectangular billiard with an inner Gaussian 
surface 

This billiard is a rectangular box with an inner Gaussian surface defined as z = 
/(r) = exp (— where a is the standard deviation and Vq is a constant 

which we will set equal to one. The classical trajectories become irregular as the 
Gaussian surface emerges in the interior of the rectangular billiard (Figure [4|. 
The trajectories were found solving the equations of motion using the fourth 
order Runge-Kutta method between successive collisions. When the particle 
collided with the frontier, then its new momentum was found and the new 
trajectory before the next collision was computed numerically. 
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(a) Eo = 0,a = 10 (b) Eo = 0,a = 0.8 (c) Eo = 0,a = 0.6 (d) Eo = 0,a = 0.4 




(e) Eo=0,a = 0.3 (f) Eo = 0,a = 0.25 (g) Eo = 10, cr = 0.25 (h) Eo = 10,a = 



Figure 4: Classical rectangular billiard with an inner Gaussian surface 

(a)-(f) The billiard in absence of an electric field, (g) and (h) The billiard is 
placed in a constant electric field. The particle has a negative charge, and it 
could not climb the surface depending if the initial kinetic energy is lower than 
energy required to move against the electric field. This case may be understood 
as a Sinai billiard with an inner soft disk. 

3.3 The phase space 

The billiards described in this document have two degrees of freedom because 
the particle moves on a two-dimensional surface. Therefore, the phase space 
is four-dimensional. Taking into account that system is Hamiltonian then the 
motion takes place on a three dimensional hypersurface of the phase space and 
we need only three variables in order to describe it. 

The Poincare section is commonly used in order to map the system and to 
obtain its dynamical information. This is equivalent to consider two variables 
which define where and how the collision occurs. These variables are the collision 
location (j) G [0,7r] and the angle a G [— f , f ] between the incident vector 
velocity projected in the plane with the normal vector of dT). The two variables 
of the mapping then define a reduced phase space. In the dunce hat billiard with 



C is a rational multiple of tt. These type of trajectories are presented in the 
reduced phase space by a set of — 1 points located in a horizontal line <p 
= constant where n' is the closed orbit period. On the other hand, if C is a 
non-rational multiple of tt, the billiard will have a quasi periodic orbit. This 
trajectory is tangent to an inner circle in the surface and it fills it generating 



0.8 



circular contour the trajectories are periodic orbits if 2/3o + eTry 1 + 
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a caustic. A quasi-periodic orbit is represented as a straight horizontal and 
densely filled line in the reduced phase space. The phase space of the dunce hat 
billiard with rectangular contour and the rectangular billiard with Gaussian 
surface were built for a sample of trajectories starting from different random 
initial conditions (see Figure [5|. When = or a ^ oo, both billiards become 
the rectangular planar billiard. The phase space of this billiard has straight 
horizontal lines typical of a non chaotic billiard. On the other hand, the phase 
space is filled by a chaotic sea when a cone or a Gaussian surface emerge in the 
rectangular billiard. 



4 Quantum Hamiltonian 



At the quantum level, the Hamiltonian of the non-planar billiards considered 
here has the form ^ 

H = -^Vl + Us + V . (32) 
The terms involved in the Hamiltonian are: 

(i) The kinetic energy on the surface. The Laplace-Beltrami operator is 



1 ^ drfd f 



1 

1 



(33) 



1 + {drfr 

with ^ plane = + + ^^0- Then, the kinetic energy may written as 

-|v|=e(r)i^„ + |.(.)a. + C(r)^, (34) 

where Hq is the free particle Hamiltonian on the plane, and K:{r) is a function 
with dimensions of wave vector defined by 

Using the equation (jsj) the radial kinetic energy on the surface |^^(^)Sr may 

be also expressed in terms of the normal Gaussian curvature as ^Kx • d^. 
Therefore, this term is null for some radially symmetric surfaces e.g. the plane 
and the cone. Finally, the non-dimensional function C,{r) is defined by C,(r) := 

(f)%(r) = l-e(r). 

(ii) The surface confining potential. The constrain of classical mechanics which 
ties the particle to the surface, in quantum mechanics may be introduced by 
a confining potential which enforces the particle to stay on the surface. This 
known as the confining approach^ and it predicts the following contribution to 
the quantum Hamiltonian [12] [13] 

U^ = ~{h-k2f (36) 

O/i 
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where ki and /c2 are the principal curvatures of the surface. For a radial sym- 
metric two dimensional surface the principal directions are along f and ^, then 
(see equation (12)) 



and 



{dr) 



_4f ^ 1# 

h dr r dr 



e(r)5. 



(37) 



(38) 



(ii) The external potential. This term includes the potential of the external 
electric field and the corresponding potential due to the boundary of the 
billiard 

V = qEofir) + V^o. (r) (39) 

where 

Jo ifr<r,((/>) 
I oo otherwise, 

and rc((/)) defines the contour of the billiard. 



(40) 



5 The quantum dunce hat billiard 



5.1 The quantum dunce hat biUiard with a circular con- 
tour 

This is probably the simplest case of study. The Gaussian curvature given 
by (|8|, is zero for this billiard as happens with the typical billiards in the plane. 
This makes a significant simplification in the Hamiltonian 



H 



"2/i 



1 



dr 



1 



Uy 



where the confining potential is 



C/s(r) 



So 



with 



R 



So • 



(41) 



(42) 



We may expect that the effect of the confining potential should not be relevant 

because it may be absorbed in — f^^^^c/) = 2^1^ least with a circular contour 

where Lz = ihdff) commutes with the Hamiltonian by virtue of the radial sym- 
metry of the billiarc|^ The analytic solution for this system may be found using 



■"^This will not be the case for an square contour where only some discrete symmetries of 
the billiard may remain depending the position of the cone center in the box. 
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the transformation defined by the equation (25). After the change of variable 
the Hamiltonian takes the form 



2/i 



1 



dr 



Td7 7 



2/if2 



(43) 



where dr = -^df-^ 9^ = V^d^ and L^j. = -^L^^. The eigenvectors and 
eigenvalues may be obtained by separation of variables. Substituting ?/^(r, (j)) = 
U{r)^{<p) in the eigenvalue equation Htlj{r^(j)) = Etlj{r^<p)^ one finds 



^{(j)) = C exp i^ifi^ 



(44) 



The angular function is ^{(j)) = C exp {im\/~^(j)) and it must fulfill with the 
condition ^{(j)) = <l>(0 + 27r), therefore fh^/^ := m G Z. As result, the angular 
solution is the same obtained in the planar circular billiard 



^{(j)) = C exp {im(j)) with m G Z. 



(45) 



The radial part of the equation will lead as usual to the Bessel differential 
equation, but with a non integer index m = m/v^, 



dr 



U{r) 



2iiEf 



2 t2 



U{r) = 0. 



Defining 



the radial solution is 



/ \2 ~ 2 zj: 

mj:{m) := m 



(46) 



(47) 



(48) 



where Prnj:{m),s (^5 = 1, 2, • • • , oo) are the zeros of the Bessel function of the first 
kind Jmj:^rn) ' "^^^ corresponding energy levels are 



^2 R'^ 

1^ f^^m-s{m),s 



2 • 



(49) 



Some eigenstates are shown in Figure [6] with a comparison with some classical 
trajectories. 

Let us consider the energy staircase function N'{E). This function counts 
the number of states below the energy E and it is defined by 



(50) 



where 0{x) is the step function. Knowing the asymptotic behaviour of the large 
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zeros of the Bessel function, one can obtain the asymptotic behaviour of Af{E) 
for large fohowing similar steps as in the planar circular billiard [14j. At the 
first order, we find that 

where Acone = ^^\J + fo ^^^^ of cone without including its base. 

The equation (51) is in agreement with the well known WeyVs formula usually 



employed for quantum billiards in the plane 



where A and P are the area and perimeter of the billiard. The confining poten- 
tial constitutes a remarkable difference between quantum non-planar billiards 
and the quantum billiards in the plane. In fact, the confining potential changes 
the spectrum and the states. Nevertheless, the agreement between the Weyl's 
formula and staircase function of the dunce hat billiard suggests that confining 
potential does not change the JV{E) at least at the first order in the asymptotic 
limit. 

Figure ([7|)-(left) shows how the energy levels are somehow "compressed" 
when the parameter fo/R is increased. Consequently, the number of energy 
levels below of a fixed energy E must increase with fo/R- This is in agreement 



with (51) 



The dunce hat billiard with a circular contour at the classical level is an inte- 
grable system. Hence, we expect that nearest neighbour energy level statistics 
will be a Poissonian in agreement with the Bohigas-Giannoni-Schmit conjec- 
ture [J. This can be verified in Figure [?]- (right). 



5.2 The quantum dunce hat biUiard with a rectangular 
contour 

By changing the circular contour of the dunce hat billiard with rectangular one, 
the rotational symmetry of the system is broken, and the angular momentum 
Lz is not conserved anymore. The system is no longer integrable. A radial 
symmetric deflector as the cone in combination with a square boundary induce 
complicated classical dynamics. Therefore, at the quantum level, the spectrum 
and eigenvectors are not simple in comparison with the ones obtained in the 
circular boundary case. We expect that the quantum dunce hat billiard in the 
square box will share some of the statistical properties in its spectrum with 
other quantum billiards with a classical chaotic counterpart. 

In general, this billiard require a numerical treatment. We implemented the 
Finite Difference Method in order to solve the eigenproblem for the Hamilto- 
nian. Some of the eigenstates are shown in Figure [S] and the nearest neighbour 
spacing distribution of the energy levels is shown in Figure |9] The cone has not 
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been placed at the center of the box in order to break the remaining discrete 
symmetries and to avoid the symmetry classification of energy levels in the level 
statistics computation. 

We have demonstrated that the staircase function of the dunce hat billiard 
with a circular box satisfies the Weyl's formula (see equation ([51])). As a result, 
we may expect that JV{E) should be linear with E for this billiard with an 
arbitrary contour, for large E. We evidence this behaviour in the numerical 
computation of M{E) (see Figure [9} (/e/^)) at least for A/'(^) < 2000. After this 
value the numerical staircase function was not linear because of the numerical 
error so we only use the first 2000 states in the level statistics computation. 

The histogram of the nearest neighbour spacing distribution for the energy 
levels below the state 2000 is shown in Figure^ {right). It fits the Gaussian 
Orthogonal Ensemble (GOE) distribution ([T]). Therefore, the change of the 
inner geometry of the two-dimensional square well from the plane to a cone 
produces the emerging of classical chaos as well as a change in the distribution 
of the energy levels of the billiard quantum counterpart. 



6 Quantum billiard with a Gaussian surface 

We may use the finite difference method in order to solve the problem for any 
surface with an arbitrary contour. For that purpose, the Hamiltonian given by 



the equation (32) is expressed as a matrix on a lattice and then it is diago- 
nalized. However, in this section, an alternative way to solve numerically the 
non-integrable problem will be presented. The idea is to use the analytic eigen- 
vectors of the planar billiard rectangular as a basis to expand the wavefunctions 
of the non-planar rectangular billiard. As an illustration, we will study the 
quantum rectangular billiard with a Gaussian surface. 



6.1 Billiard with a rectangular contour 

The eigenvectors of the rectangular plane billiard will be used as a basis ^ := 
{| -u) : = 1, 2, 3, • • • , oo} to expand the eigenfunctions of the Hamiltonian (32), 
where 

(f \u) = 4j sin (4") x) sin (4^) y) , (53) 



/A 

A = Ixly is the area of the box, and 

m{u)TT n{u)TT / / N / ^ 1 r, ^ /r .^ 
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The first kinetic term of tiie Hamiltonian (^{r)Ho expressed in ^ is 
{u\^{r)Ho\v) = I d?r{u\r)^{r){r\v) 

= e„ ^(-1)" / dV> I r-)(a./)2»{f| 



d f{u I r)r ^ exp I — \ {f\ v 



where 



(55) 



2/i 



m('u)^7r^ n('u)^7r^ 



/2 



(56) 



are the eigenvalues of the plane billiard. In (55), we have used the expansion 



1 + {drff 



(57) 



which requires | {drf)^ \< 1, i.e., a > Gmin = (^2^p^^ ^ convergence of 

the geometric series. 

Defining := ^ and using the binomial formula for r^"^ = (x+^)^^ in (55 ), 
we obtain 



{u I ar)Ho I v) 



2n 



k=0 



(58) 



where 



{an] u, v) := J x"^ exp (-a^x^^) sin (^t^^^'^^'^ sin (^/^[^'^x'^ dx' . 

The super index [s^, x^] means that the integrand has the product sin ^/^i^^^x'^ sin ^a^I^ 
and upper integration limit is l^^ with xi = x and X2 = y. All the other terms 
in the Hamiltonian may be expressed in a similar way. After a lengthly, but 
simple algebra, we find 
Confining potential terms 



{u\kUv) = 



oo , s 2(n+l) n^z 



n=0 



[ 27ra4 j 



n+2 



n+1 



/c=0 



/ . I iL i^2(n+l-/c),n^2/c,n ^ 1 ]l ] ^2{n-k) ,n2k,- 



k=0 
n+2 



^ P2(n+2-/c),n 2fe,n 



(59) 
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and 



(»i^i^)=|f:(-i)'-«(^) 



2(n+l) . n 



E 

fc=0 



I^2{n-k),n2k,n 



where the functions / must be evaluated at (a(„+i); u, v), and 

(m) ^ (m + n- 1)! 
(m-l)!n! ' 



(60) 



(61) 



Radial kinetic term 



n=0 



with 



and 



6„ := (jifj) 



2(n+l) 



(62) 
(63) 



/ I 2m -CK^r^^o i \ _ ^ \^ / ^ I [ {x) t[sc,x] t[s'^ M _i_ {y) A^^ T[sc,y] 

{U\r e X'O^ \ ~ ^ U j ^2(m-/c) + l,n^2/e,n + '^v ^2(m-/c),n^2/c+l,7 



k=0 



(64) 



where 



4!n ""'^ (c^n; '?^) := ^ exp {-anx''^) sin (^^^^x') cos dx' . (65) 

Centrifugal term 



ar)Ll 



V) = - 



4 ^ft^x ^ 



A \2nJ 

^ ' / n=l 



cxj / -r^ \ 2n n— 1 



k=0 
2 



n — 1 



(K^yA /"[^^'^^ /"[^"^'^^ /^K-(^)^ /"[^^'^^ /"^^^'^^ -k 

l^^i^ y -^2(n-k),n2k,n ^ J ^2(n-l-/c),n^2(/c+l),n^ 

-hZAi^ Ai^ -^2{n-k)-l,n2k+l,n^'^v ^2{n-k)-l,n2k,n ^ 
^ ^2(n-l-/c),n^2/c+l,n 



^2(n-l-/c) 

where the functions / must be evaluated at {ai\u^v). 
Electric potential term 



(66) 



(67) 



where the functions / must be evaluated at {a\/2;u,v). 
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The functions ll%'^\a;u^v) and ll%'^\a;u^v) can be computed by succes- 
sive differentiations with respect to a of the corresponding function with (7 = 0, 
which can be computed numericahy. Using the matrix elements previously com- 
puted, the Hamiltonian is diagonalized numerically. 

The numerical results for the rectangular billiard with the Gaussian surface 



are shown in Figure 10 We have computed 5050 energy levels and the first 
3500 of them were used to compute the nearest neighbour spacing distribution. 
The Gaussian surface is located in one of the rectangle corners. We have used a 
rectangle instead of a square in order to avoid the energy level classification by 
each symmetry. As in the case of a cone, the Gaussian surface introduces chaos 
in the classical motion which is translated at the quantum scale as a change of 
level statistics from a Poisson to a GOE distribution. Although, the histogram 
does not fit with the GOE distribution it does with the distribution function 
P2(5', pi) for mixture of chaos and regularity given by [15] 

P2(^,Pi) =/>?exp(-/>i^)erfc (^^pS^^(^2p,p ^ ^np^S^ exp (^-p,S - ^np^ 

(68) 

with p := 1 — pi. The function P2{S^pi) is a Poisson distribution for pi = 
1 and it characterizes the regular behaviour of the classical counterpart. On 
the other hand, the function P2(6', pi) is a GOE distribution ([T]) for pi = 
when the classical counterpart is fully chaotic. In general the Weyl's formula 
does not apply for this billiard because the confining potential may change the 
spectral staircase function. However, the effect of this potential is negligible for 
high energy levels. This may be appreciated at the region where the numerical 
spectral staircase function has a linear behaviour (see Figure IlO]-(c)). 



6.2 Quantum billiard immersed in a strong electric field 

If the applied electric field is strong, then it would be unlikely to find the particle 
near the top of the Gaussian surface, at least for small values of the kinetic 
energy. On the other hand, the function ({r) = 1 — ^(r) affects the particle 
motion around the top of the surface where the slope of the Gaussian function 
is not zero. Therefore, in this limit the contribution of the centrifugal term is 
negligible because the particle hardly reach the region where ({r) is important. 
The situation is similar for the K:{r) term, so under the condition of large \Eo\ 
the kinetic energy will be exclusively in S^Hq- Then, the Hamiltonian takes the 
form 

H = ^{r)Ho - ^{kr - k^f + qEof(r) + Ho. {rM)) ■ (69) 

The eigenvalues and eigenvectors can be found numerically either by using the 
finite difference method, or the method explained in the previous section. 

Two billiards will now be studied: the rectangular billiard with the Gaussian 
surface and its counterpart with a circular contour, both under the influence of 
a strong external electric field. These rectangular and circular billiards are 
analogous to the Sinai and annular billiard with a soft inner disk respectively. 
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For the circular bihiard, if the center of the Gaussian surface is located at the 
center of the circular boundary, the system has a continuous rotational symme- 
try and the angular momentum Lz is conserved. This two-dimensional system, 
having two constants of motion {H and Lz) is then integrable. The nearest 
neighbour spacing distribution of the energy levels is a Poisson distribution, a 
situation similar to the one shown in Yignie^ (right) for the dunce hat billiard 
with circular contour. If the rotational symmetry is broken by not placing the 
center of the Gaussian on the center of the circular boundary, then the system 
will be chaotic. 

The symmetry group of the rectangular billiard with the Gaussian surface is 
the dihedral group D/^. This is the same symmetry group of the Sinai billiard. 
Therefore, the rectangular billiard level statistics requires the energy levels clas- 
sification by each billiard symmetry. We may avoid this symmetry classification 
by placing the Gaussian function away from the square center in order to break 
the billiard symmetries. If a symmetry axe is hold as we show in Figure [TT]-b and 
Figure [l2]-b (this also applies for the circular billiard), then the energy levels are 
divided into two sets according to the parity of the wavefunction. The nearest 
neighbour spacing distribution of the whole spectrum is the superposition of 
two independent GOE distributions, known as the G0E2 distribution 

^ / N 1 f s^tt\ its ( s^i\\ , ( \/7rs\ , 
Pao.2{s) = -exp (--j + -exp (-^jerfc (^j . (70) 

The histograms, for both the circular and rectangular billiard in this situation, 
are shown in Figures [ll}d and[l2}c. 

On the other hand, if we break all the geometrical symmetries then the 
nearest neighbour spacing distribution of the rectangular billiard with strong 
field is a GOE distribution ([T]). This is a feature of classically chaotic systems 
with time reversal symmetry where the energy levels are likely to repel to each 
other (see Figure [TT]-e) . 

When one symmetry axe is kept, the GOE distribution can also be obtained 
by taking the energy levels of the odd or even states separately in the level 
statistics computation. Nevertheless, this process requires the computation of 
more energy levels because the parity classification enable us to use only ap- 
proximately a half of the numerically admissible energy levels computed for each 
nearest neighbour spacing distribution. 

Scars and bouncing ball states are another interesting phenomena observed 
in the standard quantum billiards. They are a manifestation of the classical 
features on the wavefunction of states in the semiclassical limit. The scarring of 
the wavefunction is a common characteristic of quantum billiards with classical 
chaotic counterpart. For the rectangular billiard with the Gaussian surface in 



a strong electric field, a wave function presenting scars is shown in Figure [13 
(d). This state corresponds to the closed unstable orbit presented in Figure 
[13]- (a). It is not a surprise to find scars in the rectangular and circular billiards 
with the Gaussian surface because of their similarity with the Sinai billiard 
and the annular billiard which also exhibit scars. However, here the scarring 



18 



of the wavefunction of these bihiards is a consequence of the external field and 
the modification of the billiard interior geometry. Taking into account that 
scarring does not appear in quantum billiards with regular classic analogue, then 
this phenomenon is another signature of chaos of the rectangular and circular 
billiards with a Gaussian surface with strong field at the quantum scale. 

Another evidence of the classical aspects on the wavefunction are the bounc- 
ing ball states (see Figure 13 and Figure 14). These states differ from the scars 
because they represent a set of classical stable trajectories on the wavefunc- 
tion. Thus a bouncing ball state represents a particle which has a well defined 
momentum but not a well defined position. 



7 Concluding remarks 

The classical chaotic behaviour of a particle in a two dimensional plane billiard 
depends of the billiard contour. However, the inner geometry of the billiard is 
also a determinant factor. We show that if we set a non planar surface into 
a billiard, then it may produce classical chaos which may be identified at the 
quantum scale, even with a regular contour (rectangular or circular). Some 
examples were studied in this paper namely: the dunce hat billiard and the 
billiard with a Gaussian surface with a rectangular or circular contours. We 
found that the Bohigas-Giannoni-Schmit conjecture remained valid for each of 
these billiards. Additionally, the billiards exhibited scarring of the wavefunction. 

For the circular billiard, if the center of the surface (cone or Gaussian) is 
placed at the center of the boundary, the system is integrable because it has 
two constant of motion, the Hamiltonian and the z-component of the angular 
momentum. We explicitly found analytically the spectrum and eigenvectors of 
the quantum dunce hat billiard with a circular contour. The effect of the surface 
was the rescaling of the spectrum by a global factor with respect to the one of a 
circular planar billiard. As a result, the nearest neighbour spacing distribution 
is the Poisson distribution. 

To confine the particle to a non-planar surface, a confining potential should 
be added at the quantum level. This confining potential affects the billiard 
energy levels. However, if the particle energy is high, then the contribution 
of the confining potential term may be dropped. As a result, although the 
Weyl's formula does not consider the quantum confining potential, it may be 
used in the asymptotic limit. This fact was analytically demonstrated with the 
quantum circular dunce hat billiard and numerically observed in the dunce hat 
and Gaussian billiards with a rectangular contour. 

The finite difference method and the expansion method were implemented in 
order to diagonalize the Hamiltonian of the quantum rectangular billiard with 
the Gaussian surface. The advantage of the second method is that a comparison 
of the numerical staircase function of both methods showed that the second one 
required the truncation of a Hamiltonian to a smaller dimension in order to 
obtain a more energy levels with an acceptable numerical error. For this aim, we 
use the fact that the slope of the spectral staircase function is just Asurfacej (47r) 
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for a large energy where the confining potential may be neglected. 

This work was supported by Facultad de Ciencias de la Universidad de los 
Andes, and ECOS NORD/COLCIENCIAS-MEN-ICETEX. 
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(a) ^ = 



: 0.1 




(c) ^ — 0.45 (d) Classical trajectory on the 

rectangular dunce hat billiard. 




(e) cr/lx = cr/ly = 0.3. (f) Classical trajectory on the 

rectangle billiard with Gaus- 
sian surface 



Figure 5: The reduced phase space. Each reduced phase space was built 
using 30 randoDi initial conditions for the first 200 collisions with the boundary. 
The phase space of the rectangular dunce hat billiard is shown in (a)-(c) for 
different values of the ratio (d) Classical trajectory in the rectangular dunce 
hat billiard setting ^ = 0.45. The phase space and a typical classical trajectory 
of the rectangular billiard with Gaussian surface setting ^ = are shown in (e) 
and (f) respectively. 
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(a) Top view of a caus- (b) Top view of the (c) State 157 

tic state 157 




Figure 6: States of the dunce hat biUiard with a circular contour. Upper 
panel Several wavefunctions of this integrable billiard exhibit the corresponding 
structure of classical dynamics. Lower panel from the left to the right. Other 
states related to the caustics. 



I 1 < state <2000 f„=R | 




Figure 7: Dunce hat billiard spectrum with a circular contour, (left) 
Changing of the billiard spectrum with the parameter fo/R- (right) Nearest 
neighbour spacing distribution of the dunce hat billiard energy levels. The dot- 
dashed, dashed and solid lines correspond to a Poisson, G0E2 and GOE distri- 
butions respectively. The nearest neighbour spacing distribution of the dunce 
hat billiard spectrum with a circular contour fits to the Poisson distribution. 
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(a) state 1 (b) state 2 (c) state 3 (d) state 4 




(e) state 5 (f) state 6 (g) state 7 (h) state 20 



Figure 8: Eigenstates of the dunce hat billiard in a square box. The 

figures from (a)-(g) are the first seven states of the non integrable dunce hat 
bihiard (top view). The symmetries of the wavefunction were broken because 
the cone was not placed at the center of the box. 



I 1 < state <3500 | 




1000 2000 3000 4000 5000 6000 
EI\Ei\ 



Figure 9: Nearest neighbour spacing distribution for the dunce hat 

billiard in a square box. (left) The numerical staircase function (solid line) 
remains linear below 2000 states, (right) Using this number of states the 
histogram of the level spacing was built. The total number of states computed 
with the finite difference method was 15625. The dot-dashed, dashed and solid 
lines correspond to a Poissonian, G0E2 and GOE distributions respectively. 
The nearest neighbour spacing distribution of the dunce hat billiard spectrum 
with a rectangular contour fits the GOE distribution. 
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(a) State 201 



I 1 < state <3500 | 




55i 6000 
^ 4000 
2000 




1000 2000 3000 4000 5000 6000 
EI\E,\ 



(b) Nearest neighbour spacing (c) spectral staircase function 
distribution 



Figure 10: Results for the rectangular billiard with a Gaussian surface. 

(a) Absolute value of the wavefunction at the state 200. (b) Nearest neighbour 
spacing distribution of the billiard energy spectrum. It fits with the distribution 
P2{S^pi) setting pi = 0.09 (dashed line), (c) The numerical staircase function 
(solid line) remains linear below 3500 states. 




(a) Eight discrete symme- (b) One discrete symmetry (c) Bihiard without symme- 
tries try axes 




(d) Nearest neighbour spacing 
distribution holding one symme- 
try axes 



(e) Nearest neighbour spacing 
distribution without symmetry 



Figure 11: Eigenstates and level statistics of the quantum rectangular 
billiard with a Gaussian surface in a strong electric field. 
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(a) Continuous rotational sym- (b) One discrete symmetry (c) Nearest neighbour spacing 
metry distribution holding one sym- 

metry axe 

Figure 12: Eigenstates and level statistics of the quantum circular bil- 
liard with a Gaussian surface with a strong field. 




(a) closed unstable orbits (b) set of classical stable tra-(c) set of classical stable trajec- 
jectories tories 





(d) scarred state 



(e) bouncing ball state 



(f) bouncing ball state 



Figure 13: Wavefunction classical features of rectangular billiard with 
Gaussian surface in a strong electric field. Upper panel Classical tra- 
jectories on the billiard. Lower panel. The corresponding scar (left) and two 
bouncing ball states. 
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(a) (b) (c) (d) 




(e) (f) (g) (h) 




Figure 14: Scars and bouncing ball states of the rectangular and cir- 
cular billiards with Gaussian surface immersed in a strong electric 
field. 
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